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1. Introduction 

The deuteron has been an object of intensive investigations as the sim- 
plest bound neutron-proton system. Throughout more than 40 years many 
methods for the description of the deuteron have been elaborated p[]-|17]. 
The main goal of any approach is the description of the interaction between 
two nucleons. It is presented either by potentials in the Schrodinger equa- 
tion or interaction kernels in the Bethe-Salpeter (BS) equation. They are 
either constructed from the presentation of the interaction as an exchange 
by various mesons (realistic potentials [IHl US]) or presented as model func- 
tions whose parameters are found from the description of observables in the 
elastic neutron-proton (np) scattering (separable ^20j-[26], phenomeno logical 
[27j models etc). There are also alternative approaches [7l[l5l[T6], relativistic 
quantum mechanics, which are based on the Hamiltonian approach. 

From the viewpoint of simplicity of performing calculations presentation 
of interaction in a separable form is the most convenient instrument [T7] . 
That is why there are separable approximations intended only to reproduce 
the behavior of the corresponding realistic potentials and used in calculations 
instead of more complicated originals (see, for example, [18] and [25]). How- 
ever, the construction of presentations of this type is a complicated problem. 
The first elaborated models [201 EI] were nonrelativistic and, therefore, they 
were of little use for the description of reactions with high-energy particles. In 
addition, there were problems with their off-shell behavior, see, for example 
[22l [23] . This behavior was adjusted in subsequent models [21]- [26] by fit- 
ting to the corresponding realistic potentials using the Ernst-Shakin-Thaler 
method [2S]. However, all these models do not contain the zero component 
of the momentum of considered particles which is necessary to construct a 
covariant model. One of the attempts of this kind is [29] where the relativis- 
tically generalized version of Graz II potential [21] was proposed. It describes 
the experimental data for the laboratory energies of the colliding neutron and 
proton TLab up to 0.5 GeV. However, its application is limited in principle 
because of nonintegrability of expressions containing the constructed form 
factors \17\ at higher TLab- The problem can be solved by using the modified 
form factors [30]. This idea was developed in [31]- [M] for the description 
of uncoupled partial-wave states in the elastic np scattering for TLab up to 
3 GeV. 

In the present paper the separable interaction kernel for the triplet partial- 
wave state ^Si-^Di is proposed. The work is a continuation of the previous 
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one [31] where the uncoupled partial-wave states with the total angular mo- 
menta J = 0, 1 were considered. Investigations of various deuteron character- 
istics are performed using the elaborated kernel. Parameters of the model are 
defined from the calculations of experimental data for phase shifts taken from 
the SAID program (http://gwdac.phys.gwu.edu) and low-energy characteris- 
tics. The off-shell behavior of the kernel is compared with the deuteron wave 
function and the Noyes-Kowalski function [35l [36] obtained for some realistic 
potentials (here Paris, CD-Bonn) which are very good at low energies. 

The paper is organized as follows. In. Sec. 2, the general Bethe-Salpeter 
formalism used for the description of the np system is considered. In Sec. 3, 
the solution of the BS equation using the separable presentation of its kernel 
is discussed. The proposed model is expounded in Sec. 4. The calculations 
of introduced parameters and obtained low-energy deuteron characteristics 
are presented in Sec. 5. The review of the results for phase shifts of the 
elastic np scattering, components of the Noyes-Kowalski function, deuteron 
wave function and the conclusions are given in Sec. 6. 

2. Bethe-Salpeter formalism 

In the relativistic field theory, elastic nucleon-nucleon (NN) scattering can 
be described by the scattering matrix T which satisfies the inhomogeneous 
Bethe-Salpeter equation [37]. In momentum space, the BS equation for the 
T matrix can be (in terms of the relative four-momenta p' and p and the 
total four-momentum P) represented as: 

T{p\p-P) = V{p\p-P) + ^^ j d'kV{p',k;P)S2{k;P)T{k,p;P), (1) 

where V{p',p; P) is the interaction kernel and 5*2 (fc; P) is the free two-particle 
Green function 

S^\k;P) = (iP-7 + fc-7-m)^'^(iP-7- A;-7-m)^'\ 

7 are the Dirac gamma-matrices. The square of the total momentum s = 
{pi + P2y and the relative momentum p = {pi — P2)/2 [p' = {p'l — P^)/'^ 
defined via the nucleon momenta pi, p2 [pi, P2] of initial [final] nucleons. 

Performing the partial- wave decomposition (see details in [13 [31]) of the 
T matrix and interaction kernel V we can rewrite the BS equation for the 
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off-shell partial-wave amplitudes: 

^a6(Po, Ip1;po, |p|;s) = K6(Po, Ip1;po, IpI; s) (2) 



+00 

I 

47]-3 

cd 



^ I dkQ /k2d|k|Kc(Po, |p'|;A;o, |k|;s) 







X 5'cd(fco, |k|; s) Tdbiko, |k|;po, |p|;s)- 

Here indices a, 6 etc denote the corresponding partial-wave state quantum 
numbers \aM) = |7r, "^^^^LjM) [3H], where S is the total spin, L is the or- 
bital angular momentum, and J is the total angular momentum with the 
projection M; relativistic quantum numbers p and vr refer to the relative- 
energy and spatial parity with respect to the change of sign of the relative 
energy and spatial vector, respectively. The two-spinor propagator Sab de- 
pends only on p-spin indices. The quantum number p defines the positive- 
(p = +) or negative-energy [p = — ) partial-wave states. In nonrelativistic 
models of nuclear-nuclear interactions only positive-energy states are consid- 
ered; therefore, p is superfluous to be shown as a quantum number. In the 
relativistic model the p = — states should be described in the general case. 
However, since only positive-energy partial states are considered in the paper 
below p is omitted. 

Calculating the T matrix we can connect the parameters of the BS kernel 
V with observables. For the description of the T matrix we use the following 
normalization condition in the on-mass-shell form for the triplet state: 

iSn f cos2eie2*^< - 1 i sin 2ei e*('^<+'^>) 
~ v^Vs - Am? V «sin2£ie^('^<+^>) cos2eie2^'^> - 1 

where m is a nucleon mass and Ei is a mixing parameter. In Eq.([3]) (5< = 
(5i=j_i, 5> = 5l=j+i and / denotes "^^^^Lj states for simphcity. Expanding 
the T matrix into a series of p-terms, according to [39j, 

pcoi5i{s) = -\ + %'' + 0{p'), (4) 

where 



P 



|p| = v/s/4-m2 = v/mTLab/2 (5) 



is the on-mass-shell momentum, one can derive low-energy parameters, the 
scattering length oq and the effective range tq. 
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The bound state of the two-particle system appears as a simple pole in 
the T matrix at s = Mj, with Md being a mass of a bound state, in our case 
it is a deutcron. Thus, the Bethe-Salpeter equation for the BS amplitude $ 
of the two-nucleon system with the total momentum J and its projection M 
has the following form: 

^'""(p; P) = J^Mp-^ P) j d'kViP, k; P)$-^^(fc; P), (6) 

The partial-wave decomposed amplitude $ can be written in the rest frame 
of the particles through the generalized spherical harmonic y and the radial 
(p part as: 

*a?(p; ^(0)) = E(^»^(p)^c)«/3 ^-(p^^ IpI; (7) 

a 

where P(o) = {Md, 0) is the total momentum of the NN system in its rest 
frame. Here Uc is the charge conjugation matrix. In the numerical calcula- 
tions instead of the amplitude it is more convenient to use the vertex function 
r which is connected with the BS amplitude: 

^jM{p;P)^S2{p;P)rjM{p;P), (8) 

so the radial parts of the BS amplitude and the vertex function g are 
connected by the relation: 

0a(po, IpI) = Yl ^<^bipo, \p\;s)gbipQ, IpI). (9) 

h 

To solve the equations for the T matrix and BS amplitude, one should 
use some assumption for the kernel V. In our case it is a separable ansatz. 

3. Separable model 

We assume that the interaction conserves parity, total angular momentum 
J and its projection, and isotopic spin. Due to the NN tensor force, the or- 
bital angular momentum L is not conserved. Moreover, the negative-energy 
two-nucleon states are omitted, which leads to the total spin S conserva- 
tion. The partial-wave-decomposed BS equation is therefore reduced to the 
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following form: 

Ti'i{Po,\p'\;po,\p\;s) = Vin{Po,\p'\;Po,\p\;s) 

+00 00 



(10) 



4vr3 



J2 dko k^c/lkl 



Vi'i"iPo, |p'h fep^ s)Ti:>i{ko, \k\;po, |p|; s) 
-E^ + lef - kl 



-00 

Supposing the separable (rank-A^) ansatz for V: 

N 

Vi'i{Po,\p'\;Po,\p\;s) = ^ Xij{s)gf\p'f^,\p'\)gf{po,\p\), 



(11) 



where the form factors Qj represent the model functions, we can obtain the 
solution of Eq. ( 10 ) in a similar separable form for the T matrix: 



N 



TniPoAp'\;poAp\;s) = J2^^A-m\poAp'\)9fiPoAp\), (12) 



where 



(13) 



hij{s) 



+ «e)2-A;2' 



(14) 



Xij{s) is a matrix of model parameters with the symmetry property: 

K,{s) = \j,is). (15) 



Using the separable ansatz (11) and performing the partial- wave decom- 
position for the interaction kernel and vertex function we can represent the 
radial part of the latter as follows: 



N 



9i{Po, IpI) = ^ij(^)9i (Po, \p\)cj{s) 



(16) 
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where / = S{D) corresponds to the ^Sii^Di) wave in the deuteron, and the 
integral equation ^ is reduced to a system of linear homogeneous equations 
for the coefficients Ci{s): 

N 

- ^ hik{s)\kj{s)cj{s) = 0. (17) 
k,j=i 

The radial part of the BS amplitude has the form (see Eq.([9 

gijpo, IpI) 

(Mrf/2-^p + ze)2-p2- 



Mpo, IpI) = TT^TTT^ — , ..A2_ 2 - (18) 



The form factors gf^ used in the separable representation ( 11 ) are obtained 
by a relativistic generalization of the initially nonrelativistic Yamaguchi-type 
functions depending on the three-dimensional squared momentum |p|. Meth- 
ods of covariant relativistic generalizations and the constructed form factors 
are discussed in the next section. 

4. Construction of the covariant separable kernel 

We consider two types of a relativistic generalization of nonrelativistic 
Yamaguchi-type form factors |10] performing the following changes: 

P' ^ V = -pI + P' (19) 

or p2 g2 : g = p - (20) 

s 

and denote the resulting modified form factors of the constructed rank-six 



interaction kernels by MY6 (19), MYQ6 (20). So high rank is necessary to 
describe simultaneously phase shifts for two waves and bound state charac- 
teristics in addition to usual low-energy parameters. The form factors have 
the following form: 

= "'If I (21) 

(Po - p^ - n)^ + "I 



JS](. ^ (Pg-P')(pc2-P0 + P')' 

g's'ip) 



[S](^^ _ (Po - P^)^(Pc3 -Po + P^^^ 
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[D] 

9i 



[S] 

9k 



ip) 



9? 



[s] 

96 



9?' = .P' 



^3 



[D] 



-(Po-P^)(Pc4-Po + P 



2\2 



9h 



9f\p) 



-{pl-v') 
{pl-v'-PlY + <' 



42) 



(22) 
(23) 



(PO - P ) (Pc6 - Po + P ) 



m-p'-piir+<mpi-p' 



/3|2)' + «: 



62 J 



The functions for the model approximation MYQ6 can be obtained from 
(21)-(23) by the change — )■ Q^. The detailed discussion of properties of 



the presented form factors can be found in |32j 

Due to the structure of the separable presentation, Eqs.(21 )-(23) 
matrix H = {hij} (14) can be written as 



the 



His) 



/ hni-s) 




\ 



h22{s) 







h23is) 











hu{s) 

/l45(s) 
/l46(s) 







/l45(s) 







/l46(s) 
/i56(s) 



(24) 



The solution of the BS equation for the radial parts of the vertex func- 
tion can be expressed through the introduced form factors with the help of 
Eqs.(|9]), (16), (17). To fix the coefficients q which are solutions of a system 
of linear homogeneous equations, we use the normalization condition for the 
^Si and ^Di states: 

Ps+pD = l, (25) 
where the pseudoprobabilities of these waves are introduced: 

The resulting coefficients q are presented in Table [T} 

The half-off-shell behavior of the interaction is controlled by the Noyes- 
Kowalski function [351 ES] in the pair rest frame (cm.) where the corre- 
sponding partial-wave decomposition of the T matrix [17] is performed: 



fi'iip; p) 





P 


;0,|p 


■,s) 


TiniO, |p 


;0,|p 





(27) 
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and compared with the Paris p[B] potential and relativistic Graz II model 
|29j . Here p (|5| is the on-shell and p is the off-shell momenta, respectively. 
It should be noted that there are three different variants of the Graz II 
separable presentation of the BS kernel differing by the D wave probability. 
We choose the variant with po = 5%. 

We also constrain the parameters of the separable interaction so that our 
model reproduces the deuteron asymptotic DjS ratio 

PD/S - 77. -7, l^oj 

where p*^ = —mEd- Finally, the deuteron magnetic moment is not in- 
cluded in the fit and is calculated as it is. 



Tabic 1: CocfRcients q. 

MY6 MYQ6 

~ci 0.0486061704 0.0704040072 

C2 0.00225058402 0.00453017913 

C3 -0.00671506834 -0.171434026 

C4 -0.019764894 0.0080399654 

C5 0.00277839389 0.00618413506 

C6 -0.304144129 0.0682392337 



5. Calculations and results 

Using the np scattering data we analyze the parameters of the constructed 
separable models. A pole in the T matrix at the mass of the bound state 
Md. 

det|r^i(. = M2)| = (29) 
is taken into account by introducing the additional parameter tjiq: 

= • (30) 

s — niQ 
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The parameter ttiq is chosen to satisfy the following condition: 

X^\s-ml) + h,,iM',)=0, (31) 

where Md = {2m — Ed) and Ed is the energy of the deuteron. 

The calculation of the parameters is performed by using Eqs.g,^ and 
expressions given to reproduce experimental values for all available data from 
the SAID program (htt p:/ /gwdac.phys.gwu.edu) for phase shifts. The low- 
energy scattering parameters are taken from jH] . 

The calculations are performed using the Wick rotation [12]. All integrals 
are calculated numerically with the technique elaborated in ^5] . 

The introduced free parameters are found from the minimization of the 

function. The asymptotic ratio po/s is also taken into account in the min- 
imization procedure. The phase shifts for S and D waves are both included 
in one minimization function: 



i=l 

+ (a^^P - a)V(Aa-P)^ + (p^;^ - p^/sf / {^Pn%f ■ (32) 

Here n is a number of experimental points. 

The effective range tq is calculated via the obtained parameters and com- 
pared with the experimental value Tq^^. 

The calculated parameters of the separable presentations MY6 and MYQ6 
are listed in Tables [2] and [3} respectively. In Table |4| the calculated low- 
energy scattering parameters and deuteron characteristics are compared with 
the corresponding experimental values and other models (Graz II [29], Paris 
[H], CD-Bonn [36] ). 

In Figsjl] and |2| the results of the phase shift calculations are compared 
with experimental data and, in addition to the afore-said theoretical mod- 
els, with the empirical SP07 SAID solution |17|. The mixing parameter is 
pictured in FigjSj As an example of comparison with a separable model we 
take here the nonrelativistic Graz II [2l] potential model (denoted by Graz 
II (NR)). The comparison with the relativistic interaction kernel Graz II [29] 
is not presented because using it the phase shifts and the mixing parameter 
cannot be calculated in the whole energy range where these observables are 
known. As it was discussed in [TTl |33], in this case when TLab exceeds some 
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limit value depending on the parameters in separable form factors it is im- 
possible to perform numerical calculations in principle, whereas our aim is to 
compare our MY6 and MYQ6 with results of other models in a wide energy 
range. Wherever it is possible to use the relativistic Graz II at high energies 
we make a comparison with it. 

The obtained components of the Noyes-Kowalski function (27) are pre- 



sented in Figs|4jj7| The 5"- and D-state wave functions 0(po,p) (18), where 
po = Md/2 — Ep [IH], are given in Figs |8] and [9| respectively, and compared 
with the corresponding relativistic and nonrelativistic Graz II models and 
the Paris potential [18j. 

As an illustration of the behavior of the elaborated separable models when 
reactions with the deuteron are considered, in Figs, 10 12, the results of cal- 
culations of cross sections for the deuteron electrodisintegration within the 
relativistic plane-wave impulse approximation are presented. We compare 
the obtained results for the differential cross section with those obtained ear- 
lier [in] using the relativistic Graz II model [22] for three different kinematic 
cases [SniET]. 
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Table 2: Parameters of the rank-six separable model with mod- 
ified ( 19 ) Yamaguchi functions MY6. 



MY6 



All 


(GeV^ 


Al2 


(GeV^ 


A 13 


(GeV^ 


Al4 


(GeV^ 


Al5 


(GeV^ 


A 16 


(GeV^ 


A22 


/ -I 7-9 

(GeV^ 


A23 


/ -I 7-9 

(GeV^ 


A24 


/ T r9 

(GeV^ 


'^25 


(GeV 


A26 


(GeV^ 


A33 


(GeV^ 


A34 


(GeV2 


A35 


(GeV2 


A36 


(GeV2 


A44 


(GeV^ 


A45 


(GeV^ 


A46 


(GeV2 


A55 


(GeV2 


A56 


(GeV2 


Aee 


(GeV2 


mo 


(GeV) 



-126.823 
-1627.106 
-78.78723 
1255.789 
1920.741 
23.25852 
-1507.037 
-202.4665 
-1211.809 
19296.73 
-19.71478 
-4.911057 
52.90785 
-557.975 
-5.781583 
2388.451 
1481.914 
23.63605 
-47615.28 
314.5085 
1.135512 
1.350753 



/3i(GeV) 

/32(GeV) 

/33(GeV) 

/34i(GeV) 

/342(GeV) 

/35(GeV) 

/36i(GeV) 

/362(GeV) 

ai(GeV) 

a2(GeV) 

as (GeV) 

0:41 (GeV) 

042 (GeV) 

«5(GeV) 

Q;6i(GeV) 

"62 (GeV) 

Pci(GeV2) 

Pc2(GeV2) 

Pc3(GeV2) 

Pc4(GeV2) 

Pc6(GeV2) 



0.1000189 
1.2089089 
0.3884728 
0.1617235 
1.0569099 
0.5975024 
0.1000189 
0.2562457 
1.7190386 
1.1342682 
0.7779747 
0.1143112 
1.9584773 

10.71719 
0.1743705 
0.3825411 
-3.947706 
-29.997902 
3.9076391 
0.5632583 

0.278038 
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Table 3: Parameters of the rank-six separable model with 



modified (201 Yamaguchi functions MYQ6. 



MYQ6 



All 


(GeV^) 


-10.574 


Al2 


(GeV^) 


-3498.482 


Al3 


(GeV^) 


2.031008 


Al4 


(GeV^) 


98.12 


Al5 


(GeV^) 


31.318 


Al6 


//^ T 7-0\ 

(GeV^) 


-30.813 


A22 


(GeV^) 


10548.08 


A23 


(GeV^) 


-76.928 


A24 


/ T 7-9\ 

(GeV^) 


-3720.413 


^^25 


(GeV j 


37427.57 


A26 


(GeV^) 


-208.697 


A33 


(GeV2) 


0.77 


A34 


(GeV2) 


-22.51814 


A35 


(GeV2) 


-84.478 


A36 


(GeV2) 


2.1435 


A44 


(GeV2) 


895.2475 


A45 


(GeV2) 


286.5565 


A46 


(GeV2) 


19.419 


A55 


(GeV2) 


-32435.14 


A56 


(GeV2) 


-705.156 


Aee 


(GeV2) 


-21.792 


mo 


(GeV) 


1.351914 



/3i(GeV) 

/32(GeV) 

/33(GeV) 

/34i(GeV) 

/342(GeV) 

/35(GeV) 

/36i(GeV) 

/362(GeV) 

ai(GeV) 

a2(GeV) 

as (GeV) 

041 (GeV) 

«42(GeV) 

as (GeV) 

a6i(GeV) 

a62(GeV) 

Pci(GeV2) 

Pc2(GeV2) 

Pc3(GeV2) 

Pc4(GeV2) 

Pc6(GeV2) 



0.100019 
0.5978559 
0.431291 
0.210503 
0.1001783 

5.5725 
0.1120488 
0.4489458 
1.468872 
1.09121 
0.7449097 

0.1 
1.945662 

4.08 
0.2058123 
0.7052189 
-4.2035 
-13.82 
9.574 
1.152 
-0.4527 
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Table 4: The low-energy scattering parameters for the triplet ^Si-^Di state and the 
deuteron characteristics. 







rot 


Pd 


Ed 


Pd/s 








(fm) 


(fm) 


(%) 


(MeV) 




(e 


/2m) 


MY6 


5.42 


1.800 


4.92 


2.2246 


0.0255 


0. 


8500^^ 


MYQ6 


5.42 


1.768 


4.65 


2.2246 


0.0259 


0. 


8535^^ 


Graz II 


5.42 


1.779 


5 


2.2254 


0.0269 


0. 


8512 


Paris 


5.43 


1.770 


5.77 


2.2249 


0.0261 


0. 


8469 


CD-Bonn 


5.4196 


1.751 


4.85 


2.2246 


0.0256 


0. 


8522 


Exp. 


5.424(4) 


1.759(5) 


4-7 


2.224644(46) 


0.0256(4)^ 


0. 


8574^= 



Calculated within the relativistic impulse approximation with the 

positive-energy partial-wave states only. 
^Ref. Hg. 
= Ref. P]. 



6. Discussion and conclusions 

The first step in construction of the interaction kernel in a separable form 
is the description of on-mass-shell characteristics, like phase shifts, low-energy 
characteristics, and the mixing parameter. In Figs jl] and |2| the results of our 
calculations of phase shifts for the ^Si and "^Di waves are presented. Figure [l] 
demonstrates that the ^5*1 partial- wave state is well described by both MY6 
and MYQ6 parametrizations for all experimental data. The Graz II (NR) 
works for TLab ^0.4 GeV. For the ^Di state MY6 and MYQ6 provide a good 
description of existing data. Graz II (NR) shows only some correspondence 
with the data for TLab ^0.4 GeV. SP07 is good for all experimental data. CD- 
Bonn was constructed for TLab ^ 350 MeV and is perfect in this region. Its 
behavior means that other models should be used for higher energies, whereas 
the interaction kernels MY6 and MYQ6 allow to perform calculations in this 
case. 

The problem of simultaneous description of phase shifts for the coupled 
partial-wave states, the mixing parameter and the low-energy and deuteron 
characteristics is worthy of a special discussion. This subject was considered 
in detail in |23] . The authors found out that an attempt to describe well the 
mixing parameter together with other observables leads to very bad results 
for the low-energy and deuteron parameters. They become too small in com- 
parison with permissible values, in particular po < 1%. It should be noted 
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that our investigation of interaction kernels of a rank ^ 6 demonstrated that 
increasing in a rank does not result in a better description of ei. Performing 
numerical investigations we discovered the same properties. In most cases 
there is no possibility to reproduce phase shifts, low-energy parameters, and 
mixing parameter simultaneously. A certain advance was achieved in the 
CD-Bonn potential model where the mixing parameter behavior looks very 
good for kinetic energies till about 350 MeV for Nijmegen group analysis [52] . 
Although there are SAID group data [17] which differ from Nijmegen group 
ones this is the best result for the moment. However, taking into account this 
discrepancy in the analysis of the experimental data for the mixing parame- 
ter we restrict ourselves to the description of all observables but ^i. At the 
same time, for completeness we present the obtained results for ei along with 
the results of all other discussed models (Fig|3|. It can be seen that both 
elaborated models do not work at all like Graz II (NR) potential which also 
does not give any agreement with the data. SP07 agrees with the experiment 
in the whole energy range. 

The half-off-shell behavior characterized by the Noyes-Kowalski function 



(27) was not fitted in any special way, it was simply calculated as it is. In 
Figs|4|j7] all components of this function are presented. The obtained MY6 
and MYQ6 functions are similar but not identical to the Graz II potential 
functions. The same can be said about diagonal components of the func- 
tion for the realistic Paris potential, whereas the difference for non-diagonal 
components is significant. In addition, MY6 modification is not identical to 
MYQ6 one. This difference should affect observables defined by the half-off- 
shell behavior of the described interaction, especially various polarizations 



As any other phenomenological model ours can describe on-shell charac- 
teristics quite easily. However, the description of the coupled ^Si-^Di channel 
is not limited only by phase shifts and low-energy observables. It is also im- 
portant to look at properties of the deuteron BS amplitude (wave function). 
Therefore, in calculations we take into account features of the wave functions 
corresponding to "^Si and ^Di parts. In Figsjs] and [9| it is shown that the 
obtained functions are similar to other discussed models in the energy region 
where their properties play a key role. 

As an example of using the interaction kernels constructed in this paper 
we calculate the differential cross sections [19] within the simplest relativistic 
plane-wave impulse approximation for various kinematic conditions [501 \5T\ 
and compare them with the corresponding relativistic Graz II model [29] cal- 
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culations. From Figs, 10 12, it can be seen that the presented results begin 
to differ only when the inffuence of the D wave increases. However, to make 
a conclusion it is necessary to take into account other effects, like the final 
state interaction (FSI) and two-body currents. This was impossible before 
because of problems with calculations, as it was discussed in [31] . Now these 
difficulties are obviated and it is planned to perform calculations of observ- 
ables with FSI using the elaborated model and the results of our previous 
work |3l] in near future. 

In conclusion, it can be said that the constructed rank-six interaction 
kernels were successfully used for the description of the on-shell and off-shell 
characteristics of the triplet ^Si-^Di partial-wave state of the np system and 
the deuteron. Good agreement with results of other models which work at 
low energies for phase shifts, low-energy parameters, deuteron wave func- 
tion was achieved. The demonstrated half-off-shell behavior is similar to the 
corresponding Graz II model. The constructed separable model of NN in- 
teraction can be used in calculations of various reactions with the deuteron, 
e.g., the deuteron photo- and electrodisintegration etc. It is also interesting 
to investigate the elastic electron-deuteron scattering process using this new 
model. However, it is a subject of a separate work. 

Additional parameters a provide integrands containing form factors of 
the separable presentation to have no poles for p component. Therefore, in 
particular, using this type of functions for form factors will make numerical 
calculations of the electrodisintegration far from the threshold possible with- 
out resorting quasipotential or nonrelativistic approximations. A comparison 
with other separable and realistic potential models allows us to demonstrate 
merits of separable kernels with a-modified form factors. For example, the 
CD-Bonn potential, which was constructed for TL^b ^ 350 MeV and works in 
this energy interval very well, cannot just be simply extrapolated at higher 
energies. As for the Graz II interaction kernel, the increase in energies of the 
considered particles is impossible in principle. On the contrary, our model 
has no these limitations. We do not describe the mixing parameter; never- 
theless, all low-energy and deuteron characteristics are reproduced well and 
the results for phase shifts cover the whole energy region and are of high 
quality. Finally, we have no insuperable obstacles in calculations and there 
is an opportunity to improve the model in future. 
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Figure 1: The model phase shifts for the "^Si and for two relativistic separable kernel cases 
MY6 and MYQ6 are compared to those of Graz II (NR) [24 , CD-Bonn \M, and of the 
empirical SP07 SAID solution gT]. 



17 



20 



O) -20 
■D 



a 

CO 



-40 



-60 



-80 



MYQ6 
MY6 

Graz II (NR) 

CD-Bonn 

SP07 



1 



Figure 2: As in Figjl] but for the "^Di wave. 
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Figure 3: As in Fig|T] but for the mixing parameter. 
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Figure 4: The Noyes-Kowalski function /oo in c.m. for two relativistic separable kernel 
cases MY6 and MYQ6 is compared to those of the separable Graz II (NR) j24] and realistic 
Paris jl8^ potentials. 
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Figure 5: As in Fig|4] but for the component /02 of the Noyes-Kowalski function. 
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Figure 6: As in Fig|4] but for the component /20 of the Noyes-Kowalski function. 
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Figure 7: As in Fig|4] but for the component /22 of the Noyes-Kowalski function. 
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Figure 8: The wave function for the '^Si partial-wave state in the deuteron rest frame for 
the MY6 and MYQ6 models in comparison with those of Graz II (NR) ^4], Graz II [29] 
and Paris [18^ . 
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Figure 9: As in Figjsj but for the ^Di partial-wave state. 
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Figure 10: Cross section for the d(e,e'p)n reaction, kinematical conditions set I [50]. n 
is the recoil momentum of the final neutron. The detailed discussion of the considered 
observable can be found in [49,1 . 
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Figure 11: Cross section for the kinematical conditions set II [SU]. See also the caption of 
FiglU 
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Figure 12: Cross section for the kinematical conditions [5T]. See also the caption of Fig 10 
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